J/ijj production and elliptic flow in relativistic heavy-ion collisions 
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Using a two-component model for charmonium production, which includes contributions from 
both the initial hard nucleon-nucleon scattering and from the regeneration in the quark-gluon 
plasma, we study the nuclear modification factor Raa and elliptic flow V2 of J/tp in relativistic 
heavy ion collisions. For the expansion dynamics of produced hot dense matter, we introduce a 
schematic fireball model with its transverse acceleration determined from the pressure gradient in- 
side the fireball and azimuthally anisotropic expansion parameterized to reproduce measured vi 
of light hadrons. We assume that light hadrons freeze out at the temperature of 120 MeV while 
charmonia at 160 MeV, similar to the kinetic and chemical freeze-out temperatures in the statistical 
model, respectively. For the properties of charmonia in the quark-gluon plasma, we use the screening 
mass between their charm and anticharm quarks and their dissociation cross sections given by the 
perturbative QCD (pQCD) in the leading order and up to the next-to-leading order, respectively. 
For the relaxation time of charm and anticharm quarks in the quark-gluon plasma, we also use 
the one calculated in the leading order of pQCD. Modeling the effect of higher-order corrections in 
pQCD by introducing multiplicative factors to the dissociation cross section of charmonia and the 
elastic scattering cross sections of charm and anticharm quarks, we find that this effect is small for 
the Raa of J/tp as they suppress the number of initially produced J/tp but enhance the number 
of regenerated ones. The higher-order corrections increase, however, the «2 of J/tp. Our results 
suggest that the V2 of J /tp can play an important role in discriminating between J/tp production 
from the initial hard collisions and from the regeneration in the quark-gluon plasma. 

PACS numbers: 



I. INTRODUCTION 

One of the signatures suggested for the formation of the 
quark-gluon plasma (QGP) in relativistic heavy ion col- 
lisions is the suppressed production of J/tp [1]. This was 
based on the idea that once the QGP is formed in heavy- 
ion collisions, the J/tp produced from initial hard scat- 
tering would be dissociated in the QGP because of the 
color screening between its charm and anticharm quark, 
reducing thus the number of produced J/tp. However, 
recent results from Lattice Quantum Chromodynamics 
(LQCD) calculations have suggested that J/tp can exist 
above the critical temperature for the phase transition 
between the QGP and the hadron gas (HG) 0,3- This 
thus not only reduces the dissociation of J/tp but also 
can lead to the production of J/tp from the recombina- 
tion of charm and anticharm quarks in QGP 4]. The 
latter effect becomes more important when the collision 
energy is high as more charm and anticharm quarks are 
produced from initial hard scattering. 

Experimentally, suppressed production of J/tp in 
heavy-ion collisions compared with that is expected from 
initial hard scattering has been observed at both Super 
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Proton Synchrotron (SPS) and Relativistic Heavy Ion 
Collider (RHIC) [^-Q. Many theoretical explanations 
have been proposed to understand these experimental ob- 
servations [9l-tl3|. Among them are the two-component 
model and the statistical model [Hj]. In the two- 
component model, some of the initially produced J/tp 
were found to survive during the collision if the J/tp dis- 
sociation temperature is sufficiently high [10]. In the sta- 
tistical model, it is assumed, however, that all initially 
produced J/tp are dissociated in QGP and the detected 
J / tp are entirely due to regeneration from charm and an- 
ticharm quarks in the QGP [l2j. Both models have been 
able to successfully reproduce the measured nuclear mod- 
ification factor Raa of J I ip as a function of the number 
of participants flClll2. 13 1 as well as of the transverse mo- 
mentum of J/tp 14j, |l5|. In the two-component model, 
the number of regenerated J/tp is smaller than that in 
the statistical model, because not all charm quarks pro- 
duced in a collision are thermalized in QGP, the finally 
produced J/tp include, however, those produced from ini- 
tial hard scattering and having survived thermal decay in 
the produced hot dense matter. Although the transverse 
momentum of thermally produced J/tp in the statistical 
model is smaller than that of J/tp produced from ini- 
tial hard scattering, it is enhanced by the transverse flow 
of J/tp. Therefore, it is difficult to discriminate the two 
models based on the number of finally produced J / tp and 
their transverse momentum distribution. 

Previously, two of us have studied in the two- 
component model the Raa of J/tp as a function of the 
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number of participants in heavy-ion collisions [13|. In 
that study, the initial local temperature of produced 
QGP was determined from the local entropy density, 
which was taken to be proportional to the linear com- 
bination of the number of participants and the number 
of binary collisions in that region. The initial fireball thus 
separates into two regions with temperatures higher and 
lower than the dissociation temperature of J/?/'- While 
the survival rate of initially produced J/ip vanishes in the 
high temperature region, its value in the low temperature 
region is finite and depends on the temperature. The dif- 
ferent values of the J/tp survival rate in the two regions 
were found to give rise to a kink in the dependence of 
the Raa of J/ip on the number of participants as seen 
in the experimental data. In the above mentioned work, 
both the screening mass between charm and anticharm 
quark and the relaxation factor of charm quarks in QGP 
were calculated in the leading order in pQCD, while the 
dissociation cross section of charmonia in QGP were cal- 
culated up to the next-to-leading order in pQCD. Ac- 
cording to LQCD, the strong coupling constant remains 
large up to around 2 T c , so below this temperature the 
QGP is nonperturbative and higher-order corrections in 
pQCD may not be negligible. The higher-order correc- 
tions are expected to enhance in QGP both the thermal 
decay rate of charmonia, which would decrease the num- 
ber of initially produced J /if), and the relaxation of charm 
quarks, which would increase the number of regenerated 
J '/if). The cancelation of these two effects is probably one 
of the reasons that both the two-component model and 
the statistical model can reproduce the measured Raa 
of J/tp as a function of the number of participants in 
heavy-ion collisions. 

In the present work, we improve our previous two- 
component model in two ways. Firstly, instead of assum- 
ing that the transverse acceleration of produced fireball 
is constant and then vanishes at a certain time as in the 
previous work, the acceleration in the present work is 
obtained from the pressure of the fireball and continues 
to be present until the fireball freezes out. Secondly, we 
take into account the azimuthally anisotropic expansion 
in noncentral collisions by using a simple parametriza- 
tion with a transverse-momentum-dependent flow veloc- 
ity to reproduce the measured «2 of light hadrons. This 
schematic fireball expansion model is then applied to 
study the Raa of J/ip as a function of both the num- 
ber of participants and the transverse momentum of J/ip 
as well as to study the V2 of J/ip- For the chemical and 
thermal freeze-out temperatures, we follow those used in 
the statistical model. In this model, the chemical freeze- 
out occurs at temperature around 160 MeV, after which 
chemical equilibrium is reached and no additional effec- 
tive chemical reactions are present and the kinetic freeze- 
out takes place later when the temperature drops to 120 
MeV. Since the kinetic freeze-out depends on the elas- 
tic scattering cross sections of hadrons and unlike light 
hadrons, the elastic scattering cross sections of J/ip in 
hadronic matter are small, we assume that the kinetic 



freeze-out temperature for J / ip is the same as the chemi- 
cal freeze-out temperature. Finally, we take into account 
the higher-order corrections by multiplying the pQCD 
cross sections with arbitrary constants to study their ef- 
fects on the Raa and «2 of J/ip. As shown later, al- 
though the fraction of J/tp produced from regeneration 
increases with increasing larger higher-order corrections, 
the resulting Raa of J/ip remains essentially unchanged. 
On the other hand, the V2 of J/ip at low transverse mo- 
mentum becomes more negative with larger higher-order 
corrections. Therefore, studying the V2 of J/ip can help 
to discriminate between the two J/ip production mecha- 
nisms of from initial hard scattering and regeneration in 
the QGP. 

This paper is organized as follows. In Sec. [Til we briefly 
describe the two-component model used in Ref. [l3| . We 
then introduce in Sec. Hill the schematic model for the 
expanding fireball produced in heavy-ion collision to ob- 
tain the time evolution of the thermal properties as well 
as the spatial and kinematic information of the fireball 
for both central and non-central collisions. In Sec. IIV1 
the Raa of J/ip as a function of the number of partic- 
ipants and the transverse momentum of J/ip are calcu- 
lated with the schematic model. Both the Cronin effect 
for the initial transverse momentum broadening and the 
leakage effect on initially produced J/ip are considered. 
Results on the elliptic flow of J/ip are given in Sec.fV] In 
Sec. lVIl higher-order corrections are introduced by multi- 
plying the pQCD results with constant factors, and their 
effects on the Raa and «2 of J/ip are investigated. Fi- 
nally, the summary and conclusion are given in Sec. I VIII 
Also, details on the the azimuthal angle dependence of 
the flow velocity are given in the Appendix. 

II. THE TWO-COMPONENT MODEL FOR J/ip 
PRODUCTION 

In the two-component model for J/ip production in rel- 
ativistic heavy ion collisions, some of initially produced 
J/ip from hard scattering can survive through the colli- 
sions and contribute to the final observed J/ip in exper- 
iments, besides those regenerated from the charm and 
anticharm quarks in the produced QGP. In this Section, 
we briefly describe the method used in Ref [l3[ to evalu- 
ate these two contributions. 



A. directly produced J/tp 

The number of initially produced J/ip in a heavy ion 
collision at impact parameter b can be determined by 
using the Glauber model, that is 

Nffc = o»%A*T AA $), (1) 

where c^,^ is the cross section for J/ip production in a 
p + p collision, which has a value of 0.774 fjh per rapidity 
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near midrapidity at collision energy y/s = 200 GeV [lj|, 
and A is the mass number of colliding nuclei. The overlap 
function T AA in the above equation is defined by 

Taa = J d 2 sT A {s)T A {b- s) (2) 

in terms of the thickness function 

/oo 
dzp A (s,z), (3) 
-oo 

where p A (s, z) is the nucleon distribution in the colliding 
nucleus. 

The survival rate of initially produced J / ip in relativis- 
tic heavy ion collisions is then given by 




where to is the time the J/ip is produced from initial hard 
scattering and T c f is the time the J/ip freezes out from 
the hot dense matter. The J/ip survival rate depends on 
its thermal decay width T(T) in the produced hot dense 
matter, 

r ( T ) = £ / 7^p^ei(fc)n. i (fc,T)af ss (fc,T). (5) 

i 

In the above, the summation is over the particle species 
i, which includes the quark and gluon in QGP as well as 
the baryon and meson in HG, with its density denoted by 
rii(k, T); v Ie i is the relative velocity between the J/ip and 
the particle; and af lss is the dissociation cross section of 
the J/ip by the particle. 

Since some directly produced J/ip are from the decay 
of directly produced excited charmonium states Xc and 
ip', the survival rate of all initially produced J/ip from 
thermal decay is the linear combination of the survival 
rates for J/ip, Xc, and ip', that is 

S th (b,s) = 0.67S^(b,s)+0.25S%(b,s) 

+0.085^(6,1), (6) 

where we have assumed that 25% and 8% of J/ip are, 
respectively, from the feed-down of Xc and ip' as usually 
adopted in the literature [13] • We note that in Ref [l3| 
the dissociate cross sections of charmonia in QGP were 
calculated in pQCD up to the next-to-leading order and 
those in HG were calculated via the factorization formula 
with the Bethe-Salpeter amplitude for the bound char- 
monium states [l8l. Il9j. 

B. regenerated J/tp 

For regenerated J/ip, their number is similarly deter- 
mined as in the statistical model, that is 

= ~f 2 {nj/^H + Br(Xc^ JW)n Xc S%_ H 

+Br(tP'^J/tP)n r S^_ H }vR. (7) 



In the above, riji^, n Xc , and are, respectively, the 
equilibrium densities of J/tp, Xc, and ip' at chemical 
freeze-out when the system has volume V. The charm 
fugacity parameter 7 takes into account the fact that the 
total charm quark number is not in chemical equilibrium. 
It is determined instead by 

where Iq and I\ are modified Bcsscl functions with no 
and rih being, respectively, the density of open and hid- 
den charms, and N^ A is the total charm quark pairs 
produced from initial hard scattering given by 

N AA = a™A*T AA {b), (9) 

where u^ 5 N is the cross section for cc pair production 
in a p + p collision, which is 63.7 fib per rapidity near 
midrapidity at y/s = 200 GeV in perturbative QCD [2(| ■ 
The factor R in Eq. (JT]) takes into account the nonequi- 
librium effect of charm quarks in QGP and is evaluated 
according to 

fl=l-exp|-y drr c (T(r))j, (10) 

where To and tqgp are, respectively, the time charm 
quark pairs are produced from initial hard scattering and 
the time QGP phase ends, and T C (T) is the thermal scat- 
tering width of charm quarks in QGP and is calculated 
in Ref using the leading-order pQCD. 

III. A SCHEMATIC MODEL FOR FIREBALL 
EXPANSION 

In our previous work [l3| , it was assumed for simplicity 
that the fireball formed in heavy-ion collisions expands 
transversely with constant acceleration of 0.1 c 2 /fm until 
certain time. Besides, the expansion is assumed to be 
azimuthally isotropic even in noncentral collisions, ne- 
glecting thus the initial elliptic shape of the collision ge- 
ometry. Since the transverse acceleration of the fireball 
is expected to be a smooth and continuous function of 
time and also to be azimuthally anisotropic in non-central 
collisions, we improve the treatment of the acceleration 
of the fireball by using the thermodynamical consider- 
ation and also introduce a simple parametrization with 
a transverse-momentum-dependent flow velocity for the 
anisotropic expansion. 

A. Central collisions 

To treat dynamically the transverse acceleration of the 
fireball formed in relativistic heavy ion collisions, we con- 
sider the pressure inside the fireball. In terms of the en- 
ergy density e, entropy density s, and temperature T of 
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the fireball as well as the chemical potential //, and num- 
ber density rii of particle species i in the fireball at a 
certain time, the pressure p is given by 

p = sT - e + ^11^ (11) 

i 

The outward force due to the pressure is then 



F = P A, 



(12) 



where A is the area of the cylindrical side of the fireball. 
Letting M be the inertia mass against the transverse ex- 
pansion of the fireball, the transverse acceleration is then 



(p-p f )A 
M 1 



(13) 



where pf is the pressure of the fireball at freeze-out, if we 
take into consideration of the fact that the acceleration 
should vanish at freeze-out, as particles undergo only free 
streaming afterwards. It is assumed in the above that 
the fireball is globally in thermal equilibrium. We have 
thus ignored the fact that in more realistic hydrodynamic 
approach only local thermal equilibrium is achieved and 
each fluid element of the fireball moves collectively. As- 
suming the transverse acceleration ar acts in the local 
frame, the transverse velocity of the boundary of the fire- 
ball in the laboratory frame at "i+l"th time step can 
then be obtained from its value at "i"th time step by 
adding the effect of acceleration relativistically, i.e., 



Ai 



1 + v l Av' 



(14) 



a l T At. 



where Av 

For the thermal quantities in Eq. (|11[) . we use the 
quasiparticle picture for the QGP, which assumes that 
strongly interacting massless partons (quarks and gluons) 
can be substituted by noninteracting massive partons to 
reproduce the energy density and pressure extracted from 
LQCD. In this approach, the pressure, energy density, 
and entropy density are given, respectively, by [2lT | 

dkMT)^~B(T) 
&i 



P(T) 


\ - 9i 

^ 6tt 2 


/ d 

Jo 
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B(T) 


e(T) 


\ - 9i 

^ 2tt 2 
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p CO 

/ d 

Jo 



-cn = E 



2^T J 



■ am l!L»*ZL. (15) 



In the above, mj(T), gi, and (Mi are, respectively, the ther- 
mal mass, degeneracy factor and chemical potential of 
parton species i with ^ obtained from the flavor conser- 
vation [22]. The parton distribution function is denoted 
by 



fi(T) 



1 



exp[CEi - mO/T] ± 1 



(16) 



with the plus and minus signs in the denominator for 
quarks and gluons, respectively, and Ei = (to 2 + fc 2 ) 1 / 2 . 
For the bag pressure B(T), it is obtained from the rela- 
tion s = dp/dT such that 



B(T) = B, + Y, [ 
.• Jt 



, drrii{T) 



d Po dm t (T) 



dT 



(17) 



where T c — 170 MeV is the critical temperature for the 
QGP to HG phase transition and Bo is the bag pressure 
at T c with its value taken to be 0.1 times the energy den- 
sity of QGP at T c to make the pressure continuous at T c . 
From the above thermal quantities, the parton thermal 
masses can be determined by fitting the QGP equation 
of state from LQCD as given in details in Ref. [1J|. For 
thermal quantities in the HG phase, they are calculated 
in the resonance gas model by including all mesons of 
masses less than 1.5 GeV and all baryons of masses less 
than 2.0 GeV, and for simplicity they are assumed to be 
noninteracting. Figure [1] shows the resulting energy den- 
sity and pressure, divided by temperature to the fourth 
power, as functions of temperature. We note that the 
pressure at the thermal freeze-out temperature 120 MeV 
is 7.75 x 10" 5 GeV 4 . 
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FIG. 1: (color online) Energy density (solid line) and pressure 
(dashed line) divided by temperature to the forth power in the 
quasiparticle model for QGP and the resonance gas model for 
HG. 



The inertia mass M in Eq. (|13|) is a free parameter 
in our model and is determined by fitting the average 
transverse momenta of light hadrons at thermal freeze 
out to the measured data. Assuming particles freeze out 
at certain proper time, their number is then 



N = 
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d 3 P 
(2tt) 3 
d 3 p 



da ■ u e~ p - u/T 



da ■ p e p ' 



pu/T 



(18) 



where a^ and are the hypersurface vector of the fire- 
ball and the current density vector of particles, respec- 
tively; and n are, respectively, the 4-velocity of the 
fluid element and the number density of particles. In ob- 
taining the last step of the above equation, we have used 
the Boltzmann distribution for evaluating the number 
density and the fact that the 4-velocity of a fluid element 
is equal to the average velocity of particles in the fluid 
element. Ignoring the transverse boost and keeping only 
the longitudinal boost, the hypersurface vector becomes 



da>* = (dT/dt,0,dT/dz)Tdr]rdrd<t> 
— (cosh 77, 0, sinh rj)Tdrjrdrd(f>, 



(19) 



where 



1, t + z 
— In 

2 t- z 



v. 



hn 1 — 
2 1 - v 



For the 4-velocity u M of the fluid element, both boosts 
are included, with the transverse boost followed by the 
longitudinal boost, leading to 

= (cosh 77 cosh p, cos <fi sinh p, sin (f> sinh p, sinh r\ cosh p) , 

where p is the transverse rapidity defined as p — tanh(uT) 
with vt being the transverse velocity. The transverse 
rapidity of the fluid element is assumed to be propor- 
tional to the radial distance from the center of the fire- 
ball, based on the result from hydrodynamics (23|. For 
the 4-momentum p^ of a particle, it is 

p M = (wit cosh y, pt cos ip, pr sin ip, ttit sinh y) , 



where y and ttit — y m 2 + p T are, respectively, the ra- 
pidity and the transverse energy. In general, three dif- 
ferent azimuthal angles are required: one for the hyper- 
surface vector a^; one for the 4-velocity of the fluid 
element; and the other for the momentum vector p^. In 
central collisions, the azimuthal angle for a^ is, however, 
the same as that for u^. This is also the case in non- 
central collisions if only the radial flow is considered and 
this is assumed in present study. 

Substituting Eqs. (19]) into Eq. ([18]) leads to [H 



dN 
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TITlT 



2 7T 



dydp\ 2(2tt) 3 J q 

dAx / dtp e 



•ip I da-p e~ p - u/T 

Pt sinh p cos((p—ip) /T 



2(2tt) 3 

poo 

x / drj coshiij - y) e- mTCOshpcosh ^-yV T 



TTTLT 



dA T I a 



Pt sinh p 



T 



rriT cosh p 



T 



,(20) 



where At is the transverse area, and I n and K n are modi- 
fied Bessel functions. From Eq. (TT8"|) . the averaged trans- 
verse momentum of particles is given by 



{pt) 



J dp T PT(dN/dydp T ) 
j dp T (dN/dydp T ) 



(21) 



We find that with the inertia mass M = 570 GeV, the 
transverse velocity of fireball at freeze-out is 0.71, and 
the average transverse momenta of pions, kaons and pro- 
tons are 462, 668 and 953 MeV/c, respectively, which are 
comparable to measured values of 451±33, 670±78 and 
949±85 MeV/c for 7r + , K + and proton, respectively, in 
Au+Au collisions as ^/s = 200 GeV and of 0-5% central- 
ity [25J]. We note that the measured average transverse 
momenta of 7r~ , K~ and antiproton have similar values 
as those for their antiparticlcs. 
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FIG. 2: (color online) Transverse velocity (a) and temperature 
(b) of fireball in central collision using transverse acceleration 
from the thermodynamical approach (solid line) or constant 
acceleration (dotted line) as functions of time. 

In Fig.[2ja), we show the transverse velocity of the fire- 
ball obtained with the acceleration using the above ther- 
modynamical approach as a function of time and compare 
it with that using the constant acceleration as in our pre- 
vious study. The transverse velocity is seen to increase 
rapidly at first and then slowly in the thermodynami- 
cal approach, whereas it increases linearly with time and 
stays constant after a certain time in our previous study 
with constant acceleration. For the temperature of the 
fireball, shown in Fig.[5Jb), it evolves in time similarly in 
the two cases in spite of their different evolutions in the 
transverse velocity. In the thermodynamic approach, the 
QGP phase terminates at t — 5.4 fm/c, the mixed phase 
at t = 6.9 fm/c, and the HG phase reaches the freeze- 
out temperature at t = 18.7 fm/c. In the case of con- 
stant acceleration, the corresponding times are 5.6 fm/c, 
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7.0 fm/c, and 18.3 fm/c, respectively, and are thus sim- 
ilar. Therefore, the modified acceleration introduced in 
the present study is not expected to change our previous 
results on the Raa of J/ip as a function of the number 
of participants [13J . 



B. Non-central collisions 

In noncentral collisions, the length of the overlapping 
region of two colliding nuclei is shorter in the direction of 
impact parameter than in its perpendicular direction on 
the transverse plane, which are designated, respectively, 
as the re-direction and the y-direction. The difference in 
the overlapping lengths along these two directions leads 
to difference in the pressure gradient in these two di- 
rections, and as a result fluid velocities along these two 
directions are different as well, with the fluid velocity in 
the a;-direction higher than the one in the y-direction. To 
include such an anisotropic expansion of the fireball, we 
add a term proportional to the eccentricity of the fireball 
to the acceleration [Eq. (|13l) ] as follows 



where 



R{<f>) 



Rl 



-1/2 



(25) 



is the radial distance from the center of the fireball to 
its surface along the azimuthal angle tj>. Details in the 
derivation of the above result are given in the Appendix. 
For the transverse acceleration a? in Eq. (|22p. it is ob- 
tained in the same way as in central collisions. The con- 
stant z in Eq. (|22l) is determined by requiring that the 
final vi of light hadrons, given by the expectation value 
of cos(2</?) as 



v 2 {pt) 



J dip cos(2<p) (dN/dyd 2 p T ) 



J d<p(dN/dyd 2 p T ) 
J dA T cos(2</>) I 2 {p T smhp/T)Kx(m T coshp/T) 



/ dA T I (p T smhp/T)K x {m T coshp/T) 



(26) 



reproduce the measured data. 



a x = a T {\ + ze) 

a y = or(l-*e), (22) 
where the eccentricity e is defined as 

Ry Rx 
Ry + Rx ' 

and z is a fitting parameter to take into account the dif- 
ferent accelerations in and perpendicular to the reaction 
plane in non-central collisions. In the above, R x and R y 
are, respectively, the half length of the fireball in the in- 
direction and in the y-direction on the transverse plane. 
If the shape of the overlapping region is a perfect ellipse, 
they then correspond to the semiminor and semimajor 
of the ellipse, respectively. For simplicity, we transform 
the initial shape of the fireball into a ellipse by keeping 
the same transverse area and the same R x to R v ratio. 
The velocities in the x-direction and in the y-direction 
are then, respectively, 



Av x 



v x Av x 
f Av y 



1 + VyAVy ' 



where Av x 



''At and Av v 



(23) 



Assuming that the 



fireball keeps the elliptic shape on transverse plane dur- 
ing expansion, the fluid velocity, which is assumed to be 
along the radial direction, as a function of the azimuthal 
angle 4> is then 



v{4>) = R 3 ((f)) 



Rl 



v x sin 



(24) 
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FIG. 3: (color online) Time evolution of temperature (a), v x 
and v y (b), R x and R y (c) for Au+Au collisions at t/s — 200 
GeV and impact parameter b —9 fm. 

In Fig. [31 we show the time evolution of the tempera- 
ture of the fireball and of v x , v y , R x and R y for Au+Au 
collisions at y/s = 200 GeV and impact parameter b — 
9 fm. These results are obtained with the anisotropic 
constant z — 2.2 in Eq. (j2"2"j) and the inertia mass 
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M = 108 GeV in Eq. {T3j|. It is seen that the QGP 
phase ends at t = 2.8 fm/c, the mixed phase lasts until 
t = 3.8 fm/c, and the freeze-out temperature is reached 
at t = 11.1 fm/c. The flow velocities v x and v y are 0.51 c 
and 0.27 c at the phase transition to HG and reach 0.76 
c and 0.57 c at freeze-out, respectively. We note that the 
length of R x becomes slightly longer than that of R y at 
freeze-out. 
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FIG. 4: (color online) Elliptic flow v% of pions (solid line), 
kaons (dashed line), and protons (dotted line) for Au+Au 
collisions at y/s = 200 GeV and impact parameter b = 9 fm 
before the high pr correction (a) and after the correction (b) . 
Experimental data from collisions at 20-60% centrality [2y] 
are shown by rectangles, circles, and triangles for pions, kaons, 
and protons, respectively. 

The elliptic flow v-i of pions (solid line) , kaons (dashed 
line) and protons (dotted line) at b — 9 fm, obtained 
from Eq. (|2"6")l and using the results of Fig. [31 are shown 
in Fig. Eta). Except in the low px region, they are seen 
to deviate from the experimental data from collisions at 
20-60% centrality [26| . shown by rectangles, circles, and 
triangles for pions, kaons, and protons, respectively. The 
reason for this deviation is because high-p^ particles are 
not likely to be completely thermalized as assumed in 
our fireball model. To take into the non-equilibrium ef- 
fect, we multiply exp[— C(pr/n)] to Aw = (y x — v y )/2, 
where C = 1.14 GeV -1 is a fitting parameter and n is the 
number of constituent quarks in a hadron, as described in 
the Appendix. Figure Etb) shows the resulting v-i of light 
hadrons after this high px correction, and they are seen 



to reproduce very well the experimental data at all px- 
Since the high px correction factor modifies the Boltz- 
mann distribution in Eq. (|18p . it can be considered as 
an effective way to include the viscous effect on the fire- 
ball expansion. In the quark-counting picture for hadron 
scattering, the ratio of the baryon to the meson scatter- 
ing cross sections is 3/2. This then leads to a smaller 
viscous effect on baryons than mesons, which is consis- 
tent with the quark number dependence in the above high 
Pt correction factor. 




400 



FIG. 5: (color online) Participant number dependence of aver- 
age transverse momenta of ix + (solid line), K + (dashed line), 
and proton (dotted line) from the fireball model with cor- 
responding experimental data shown by the rectangle, circle 
and triangle. 

The average transverse momenta of pions, kaons, and 
protons in the above case of b — 9 fm are 433, 619 and 
874 MeV/c, respectively, which are comparable to mea- 
sured values of 436±35, 655±54, and 901±88 MeV/c in 
collisions at 30-40% centrality as well as 426±35, 636±54 
and 868±88 MeV/c in collisions at 40-50% centrality for 
7r + , K + and proton, respectively [2^|. For other val- 
ues of the impact parameter, it is found that the inertia 
mass parameterized as M = 0.36Np£% with N pai t being 
the number of participants in a collision, works well as 
shown in Fig. [5j where the participant number depen- 
dence of the average transverse momenta of ir + , K + and 
protons obtained with the parameterized inertia mass are 
compared with the experimental data. 



IV. NUCLEAR MODIFICATION FACTOR OF 

The nuclear modification factor Raa of J/ip is a mea- 
sure of its production in heavy ion collisions relative to 
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that in p+p collisions multiplied by the number of binary 
collisions. Its value is less than one if there is a suppres- 
sion and greater than one if there is an enhancement. 
Besides effects due to the produced hot dense matter, 
the Raa of J/ip is also affected by two other important 
effects, namely the Cronin effect and the leakage effect. 



A. the Cronin effect 

Charmonium is produced mainly through the gluon- 
gluon fusion in high-energy nuclear collisions. If a gluon 
in one nucleus is scattered by the nucleon in the other 
colliding nucleus before it produces the charmonium, the 
gluon obtains an additional transverse momentum, and 
the charmonium produced from the gluons then also ac- 
quires this additional transverse momentum. As a result, 
the distribution of the transverse momentum of J / ip be- 
comes different from that in p+p collisions, and this is 
called the Cronin effect. Taking into account of this ef- 
fect, the average transverse momentum of J/ip in A+A 
collisions is modified as following: 



(Pt)aa 



(Pt)pp 



(Sp 2 T ) 



A 



yN 



(L 9 n), 



(27) 



where (p^)aa an d (j>t)pp are average squared transverse 
momenta of charmonium in an A+A collision and in a 
p+p collision, respectively; (5pj) is the average of added 
squared transverse momentum by one gluon-nucleon col- 
lision, X g N is the mean free path of a gluon in uniform 
nuclear matter, (L 9 m) is the average path length of the 
gluon in uniform nuclear matter before it produces a 
charmonium. The prefactor (Sp?n) /A 9 jv is fitted to the 
experimental data of p+A or d+A collisions and (£ ff jv) 
is calculated with the Glauber model according to 



(L gN ) - 



J d 2 bd 2 sdzAdzBL(b, s, za, zs)K(b, s, za, %b) 
J cPbd 2 sdzAdzBK(b, s, za, Zb) 



where 

- f ZA PA (s,z) 
L(b,s,z A ,z B ) = / dz 

J— oo Po 



dz 



(28) 



p B {b-s,z) 

Po 



K(b,s, za, z B ) = Pa{s, z A )pB{b - s, z B ) 

dzpA (s, z) + / dzps{b — s,z) 

1 ZA — OO 



X CXp - Oabs 



L(b,s,ZA,ZB) is the sum of the path length in nucleus 
A of a gluon in nucleus B and the path length in nu- 
cleus B of a gluon in nucleus A, which correspond to the 
first term and the second term on the right hand side of 
L{b, s, za, zb), respectively. The probability that a char- 
monium is produced at the position za from the center of 
nucleus A and zb from the center of nucleus B is given by 
K (b, s, za, zb)- The exponential factor in K(b, s, za, Zb) 
is due to the effect of nuclear absorption or nuclear de- 
struction as produced J / ip can be absorbed or destroyed 
by nucleons with the absorption cross section a a bs before 
it gets out of the colliding nuclei. Notice that the integra- 
tion range for the Cronin effect is from — oo(oo) to za{zb) 
and for nuclear absorption from za(zb) to oo(— oo) in 
nucleus A(B), as the Cronin effect is present before the 
production of the primordial charmonium and the nu- 
clear absorption happens after the creation of the J/ip- 
In the present work, L 9 n is calculated at each point on 
the transverse plane as following: 



L gN {b,s) 



J dz A dz B L(b, s, z A , z B )K(b, s, z A , z B ) 
J dz A dz B K(b,s,ZA,ZB) 



(29) 



The Cronin effect thus depends on the transverse position 
at which a J/ip is produced. 

The transverse momentum distribution of J/ip in A+A 
collisions is related to that in p+p collisions by 



fA+A{PT,b,S) 



1 



d 2 p' T exp 



p'j 



5p 2 T {b,s) 



x fp+p(\pi 



P' T \ 



(30) 



where 



8p%(b,s) 



( S Pt) T ( t * 



A 



yN 



and fp+pipr) is the distribution function of J/ip trans- 
verse momentum in p+p collisions. For center of mass 
energy at y/sW^ = 200 GeV, f p+p {pr) = a(l +p^//3)~ 6 
with (3 = 4.1 GeV 2 0, and we use {5p 2 T )/\ gN = 
0.1 GeV 2 /fm and a a b s = 1.5 mb for Au+Au collisions 

a 



B. the leakage effect 



In the above, b and s are the impact parameter and trans- 
verse vector in configuration space, respectively; za and 
zb are longitudinal positions where J/ip is produced from 
the center of nucleus A and from the center of nucleus 
B, respectively; pa and ps are the distribution functions 
of nucleons in nucleus A and in nucleus B, respectively, 
and they are normalized to their mass numbers; p Q is 
the saturation density of nuclear matter. The function 



In our previous work, it was assumed that charmonia 
produced inside the fireball could not get out and char- 
monia produced outside could not get into the fireball. 
In realistic situation, it is possible for charmonia of high 
Pt to escape from the fireball, leading to the so-called 
leakage effect. To take into account this effect, we con- 
sider the thermal decay width of a charmonium inside 
the fireball, i.e., its transverse position (x and y) satisfies 
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[x/Rxit)] 2 + [y/Ry(t)} 2 < 1, 

s—y f d 3 k 
T(x,y,t) = Y,9i] J2^ v * 



.x{k) ni (k,T)af ss {k,T) 



(31) 



where i denotes the particle species that dissociates the 
charmonium with gi and rii being its degeneracy factor 
and number density, respectively; v rc \ is the relative ve- 
locity between the charmonium and the particle; af lss 
is the dissociation cross section of charmonium by the 
particle. The charmonium thermal decay width vanishes 
once it moves out of the fireball when its transverse po- 
sition satisfies [x/R x (t)] 2 + [y/R y {t)\ 2 > 1. Furthermore, 
we set the charmonium decay width to infinite if it is 
produced in a region with temperatures higher than its 
dissociation temperature. 
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FIG. 6: (color online) Raa of J/ip as a function of the num- 
ber of participants including both the Cronin and the leakage 
effect. The dotted line is the Raa for constant acceleration 
and the solid line is that for acceleration determined thermo- 
dynamically, which are, however, almost indistinguishable. 

Results obtained with the two-component model de- 
scribed in Sec.|TT]for the Raa of J/ip as a function of the 
number of participants including both the Cronin and 
the leakage effect are shown in Fig. [BJ The dotted line is 
the Raa with constant acceleration of the fireball and the 
solid line is that with the fireball acceleration determined 
thermodynamically. In the latter case, the Cronin effect 
and the leakage effect are included. As shown in Fig. 
Hfb), the lifetime of QGP is slightly longer in the former 
case than in the latter case, and this leads to a somewhat 
smaller Raa in the former case. Since the Cronin effect 
only changes the transverse momentum of J /if), not the 
number of J/ip, it does not affect the Raa of J/ip- The 



FIG. 7: (color online) Nuclear modification factor Raa of 
J/ip as a function of transverse momentum in collisions of 0- 
20% (a), 20-40% (b), 40-60% (c) and 60-92% (d) centralities. 
Dashed and dotted lines are, respectively, the Raa from ini- 
tially produced J /if) and regenerated J/ip, and the solid line 
is the sum of the two. 



leakage effect, on the other hand, may increase Raa or 
decrease Raa, depending on whether J/ip initially pro- 
duced outside the fireball can move into the fireball and 
on whether J/ip initially produced inside the fireball can 
escape to the outside. The net leakage effect on Raa is, 
however, found to be small, and the difference between 
the results for constant acceleration and those for acceler- 
ation determined thermodynamically is thus hardly seen 
in Fig. [5J 

Figure shows the Raa of J/ip as a function of trans- 
verse momentum obtained with the fireball acceleration 
calculated thermodynamically. It is seen that the Raa 
of initially produced J/ip, given by the dashed line, in- 
creases with transverse momentum because of the Cronin 
effect. On the contrary, the Raa of regenerated J/ip de- 
creases with transverse momentum at high px as shown 
by the dotted line. We note that the Raa of J/ip is 
slightly enhanced and suppressed at high and low px, 
respectively. 

In Figure [51 we show the participant number depen- 
dence of the average squared transverse momentum of 
J / ip in Au+ Au collisions from the two-component model. 
Our results reproduce the experimental data for periph- 
eral and mid-central collisions. For central collisions, the 
two-component model overestimates the measured aver- 
age squared transverse momentum of J/ip. One possible 
reason for this is that we may have overestimated the 
diffusion of charm quarks in the fireball. The number 
of charm quarks produced in heavy-ion collision is pro- 
portional to the number of binary collisions of nucleons, 
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FIG. 8: (color online) Average squared transverse momen- 
tum of J/tp in Au+Au collisions at ^/snn = 200 GeV from 
the two-component model (solid line) and experimental data 
(filled squares). 



which is more concentrated in the central region of the 
fireball. The number of light quarks is, on the other hand, 
proportional to the number of participants. Since the 
relative concentration of the number of binary collisions 
to the number of participants is larger in central colli- 
sions than in semi-central or peripheral collisions, charm 
quarks are initially more concentrated in the central re- 
gion than light quarks. If the diffusion of charm quarks is 
slow, they would remain near the central region even at 
hadronization and their transverse velocity would then 
be smaller than those of light quarks. 



V. ELLIPTIC FLOW OF J/tp 

The elliptic flow V2 of J/tp in the two-component model 
can be different from that of light hadrons as a result of 
following two effects. Firstly, J/tp's are produced from 
two different mechanisms: the initial hard scattering and 
the regeneration, while light hadrons are produced only 
through regeneration. Secondly, the freeze-out tempera- 
ture for J/tp could be different from that of light hadrons. 
Since initially produced J/tp's have transverse momenta 
that are azimuthally isotropic, they are expected to have 
vanishing V2 ■ However, the anisotropic expansion of the 
fireball could give rise to nonzero V2 for J/ip. In non- 
central collisions, the size of the fireball is shorter in the 
x-direction than in the y-direction. As a result, the sur- 
vival rate of J / ip moving in the x-direction will be higher 
than that moving in the y-direction. In this case, the 
v% of J/tp is positive, because the number of J/tp mov- 



ing in the x-direction is more than that moving in the 
y-direction. On the contrary, the of regenerated J/tp 
is negative at low px for the following reason. Assuming 
particles are locally thermalized, they then have a ther- 
mal distribution in the local frame. Since the thermal 
distribution is peaked at p 2 = 2T(T+ \JT 2 + m 2 ), where 
T and m are the temperature and mass of the particle, 
respectively, more massive particles thus peak at higher 
momenta. If the thermalized matter is boosted with a 
certain velocity, the peak is shifted to higher momen- 
tum in the laboratory frame. In noncentral collisions, 
the difference in the length of the overlapping region in 
the x-direction and the y-direction leads to a difference in 
the pressure gradient and as a result a difference in fluid 
velocity along these two directions. Because the length 
of the overlapping region is shorter in the x-direction, 
the fluid velocity is higher and the peak of the momen- 
tum distribution is more shifted in that direction. As 
the peak shifts to higher momentum, the number density 
at momentum lower than the peak momentum becomes 
smaller than the number density in the unshifted or less 
shifted momentum distribution, while the number den- 
sity at momentum higher than the peak becomes larger. 
As a result, the number density of J/tp moving in the 
x-direction is smaller than that of J/tp moving in the y- 
direction at small pr, and this becomes opposite at high 
Pt, leading thus to a negative V2- 
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FIG. 9: (color online) Elliptic flow V2 of J/tp for collisions at 
b — 9 fm in the two-component model. The dashed, dotted 
and solid lines are, respectively, the V2 of initially produced, 
regenerated and total J/tp. 

Figure [5] shows the results from the two-component 
model for the V2 of initially produced, regenerated and 
total J/tp in Au+Au collisions at ^/snn = 200 GeV and 
b = 9 fm. These results are obtained with the high-p^ 
correction to the flow velocity as in the case for light 
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hadrons. It is seen that the «2 of initially produced J/tp 
is essentially zero, whereas the of regenerated J/tp 
has a minimum value of —0.028 at pt = 1.1 GeV. The 
V2 of total J/ip is, however, slightly negative, because 
only 12% of J/tp are from regeneration for collisions at 
b = 9 fm in the two-component model. We note that the 
V2 of J/tp has been recently measured by the PHENIX 
collaboration for Au+Au collisions at 20-60% centralities 
[27l HI] and it has a value of -0.15±0.12 between p T = 1 
and 2 GeV, which is much much more negative than our 
value. 



VI. EFFECT OF HIGHER-ORDER PQCD 
CORRECTIONS 

In our two-component model, the screening mass be- 
tween charm quarks and the thermal decay widths of 
charmonia as well as the relaxation factor of charm 
quarks in QGP are calculated in pQCD. Specifically, the 
dissociation cross section of charmonium used for de- 
termining its thermal decay width is calculated up to 
the next-to-leading order, and the elastic cross section of 
charm quarks for calculating the relaxation factor, and 
the screening mass are calculated in the leading order. 
Since it is known that pQCD is not reliable up to around 
2 T c , corrections due to higher orders are expected to 
be important for the temperature regime considered in 
the present study. The higher-order correction for the 
dissociation cross sections of charmonia is expected to 
suppress the survival rate of charmonia in hot dense nu- 
clear matter while the correction for the relaxation fac- 
tor of charm quarks would enhance the regeneration of 
charmonia. Therefore, there could be a balance between 
the dissociation and the regeneration of charmonia in 
calculations including higher-order corrections. To in- 
clude higher-order corrections, we simply introduce in 
the present study some multiplication factors to the dis- 
sociation cross sections of charmonia and to the elastic 
cross section of charm quarks as following: 

a J/-ip+q(g)^c+c+X = ^ <J J/ip+q(g)^c+c+X 

a 'c+q{g)^c+q(g) = B(T c+q(g)^c+q(g), (32) 

where A and B are arbitrary constants. 

In Fig. [TU1 the Raa of J / ip is shown as a function of 
the number of participants. Solid, dashed, dotted, and 
dash-dotted lines correspond, respectively, to results ob- 
tained in the two-component model with (A,B) = (1,1), 
(2,2), (5,4), and (10,6.5) for the multiplicative factors 
used to include high-order effects in pQCD. These re- 
sults show that higher-order corrections do not change 
much the value of the Raa of J/tp even with a large cor- 
rection factor of ten. This is so because of the balance 
between the dissociation and the regeneration of J/tp as 
mentioned previously, with the suppression of J/tp caused 
by enhanced dissociation compensated by the enhance- 
ment of regenerated J/tp. However, as the correction 
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FIG. 10: (color online) Nuclear modification factor Raa of 
J/tp as a function of the number of participants with different 
values of the multiplicative parameteres (A,B)=(1,1), (2,2), 
(5,4) and (10,6.5) for including the effect of higher-order cor- 
rections in pQCD. Corresponding results are shown by solid, 
dashed, dotted, and dot-dashed lines, respectively. 



factor becomes large, the kink in the Raa curve around 
Apart = 170 becomes less noticeable. To understand this, 
we note that the kink in the curve is caused by the dis- 
continuity in the survival rate of J/tp in the region of 
the fireball where the temperature is higher than the dis- 
sociation temperature and in the region where the tem- 
perature is lower than the dissociation temperature. As 
higher-order corrections become large, the thermal decay 
width of J/tp increases and the survival rate of J/tp in 
the low temperature region becomes comparable to the 
survival rate in the high temperature region. As a re- 
sult, the kink disappears when higher-order corrections 
are large. 

Figure ITTI shows the fraction of regenerated J/tp among 
total J/tp in heavy-ion collisions in the two-component 
model. Without higher-order corrections (solid line), the 
regenerated J/tp is about 31% of total J/tp produced 
in central collisions. Including higher-order corrections, 
the percentage increases to 52, 81 and 96%, respectively, 
for multiplicative factors (A,B)=(2,2), (5,4) and (10,6.5), 
shown by dashed, dotted, and dot-dashed lines, respec- 
tively. We note that the faction of regenerated J/tp is 
smaller in semi-central and peripheral collisions than in 
central collisions. 

Figure [T2] shows the Raa of J/tp as a function of the 
transverse momentum in collisions of 0-20% (a), 20-40% 
(b), 40-60% (c) and 60-92% (d) centralities. The lines in 
each figure correspond, respectively, to multiplicative fac- 
tors (A,B)=(1,1), (2,2), (5,4) and (10,6.5) for the higher- 
order corrections. As we can see, Raa does not change 
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FIG. 11: (color online) Fraction of regenerated J/ip among 
total J/ip in heavy-ion collisions. Solid, dashed, dotted, and 
dot-dashed lines correspond, respectively, to multiplicative 
factors (A,B)=(1,1), (2,2), (5,4) and (10,6.5) for higher-order 
corrections. 
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FIG. 12: (color online) Raa of J/ij) as a function of transverse 
momentum in collisions of 0-20% (a), 20-40% (b), 40-60% (c) 
and 60-92% (d) centralities. Solid, dashed, dotted, and dot- 
dashed lines correspond, respectively, to multiplicative factors 
(A,B) = (1,1), (2,2), (5,4) and (10,6.5) for the higher-order cor- 
rections. 



much even with large higher-order corrections as a re- 
sult of the large transverse flow of regenerated J/ip. As 
shown in Fig. the px dependence of the Raa of ini- 
tially produced J/ip is similar to that of regenerated J/ip 
except at high where the Raa of initially produced 
J/tp increases while that of regenerated J/ip decreases. 
As a result, the Raa at high px decreases when higher- 
order corrections are included, and this seems to contra- 
dict with the pt dependence of measured Raa, which is 
an increasing function at high pt- Although higher-order 
corrections are small at high pt, they are relatively large 
at low pt, and this seems to reproduce the observed de- 
creasing Raa with increasing pt, when pt is low, in pe- 
ripheral collisions as shown in Fig. I12f d). This is mainly 
due to the fact that in peripheral collisions, the time 
interval between the thermalization of fireball and the 
frcczc-out of J/ip is short, so the transverse flow is not 
sufficiently developed to give a large transverse flow ve- 
locity to the regenerated J/ip. 

In Fig. 1131 we show the average squared transverse mo- 
mentum of J/ip before and after including higher-order 
corrections. It is seen that higher-order corrections have 
very little effect on the results. 

Figure [14] shows the V2 of J/ip in collisions at b = 9 
fm without and with higher-order corrections. Different 
from the Raa of J/ip, the change in V2 due to higher- 
order corrections is significant. For increasing large cor- 
rection factors listed in the above, the fraction of regen- 
erated J/ip among total J/ip is 12%, 23%, 42% and 63% 
respectively. This can be easily understood from the fact 
that the vi of initially produced J/ip is very small. In 



terms of the v% of initially produced J/ip and that of 
regenerated J/ip, the vi of total J/ip is given by 

Jcos(2tp)dN f cos(2(p)dN R + fcos(2<p)dNi 
V2 ~ JdN ~ J dN R + J dNi 

/ cos{2<p)dN R J cos(2y)dN R N R 
~ J dN R + J dNj JdN R N ' [ ' 

where N, Nj and N R are, respectively, the total num- 
ber of J/ip, the number of initially produced J/ip and 
that of regenerated J/ip. Equation (f33|) shows that the 
v-2 of total J/ip in heavy-ion collisions can be approx- 
imated by the v% of regenerated J/ip multiplied by its 
fraction among total J/ip. The double-dot-dashed line 
in Fig. [14] is the V2 of J/ip if all J/ip are produced from 
regeneration. It is worthwhile to point out that even in 
this extreme case, the V2 of J/ip still only reaches the 
upper error bar of experimental data, which has a value 
of -0.03 between p T = 1 and 2 GeV [13, H|. If the ex- 
perimental data on V2 of J/ip is correct, it then means 
that the fraction of regenerated J/ip is not small even in 
semi-central collisions. 



VII. SUMMARY AND CONCLUSION 

We have investigated in the two-component model the 
Raa of J/ip as a function of transverse momentum and 
as a function of the number of participants as well as 
the V2 of J/ip in relativistic heavy ion collisions. Firstly, 
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FIG. 13: (color online) Average squared transverse momen- 
tum of J/ip in Au+Au collision. Solid, dashed, dotted, and 
dot-dashed lines correspond, respectively, to multiplicative 
factors (A,B) = (1,1), (2,2), (5,4) and (10,6.5) for the higher- 
order corrections. 

a schematic model was constructed to describe the ex- 
pansion of the fireball formed in these collisions. The 
transverse acceleration of the fireball was assumed to be 
proportional to the outgoing force given by the product 
of the pressure of the hot dense matter inside the fireball 
and the surface area of the fireball, and inversely pro- 
portional to an inertia mass, which was taken to be a 
free parameter to fit the experimental data on the mea- 
sured transverse momentum spectra of light hadrons. 
The anisotropic expansion in non-central collisions was 
modeled by simply adding or subtracting an additional 
term proportional to the eccentricity of the fireball to the 
acceleration in x-direction or y-direction. The coefficient 
of the additional term was determined such that the V2 of 
light hadrons was reproduced after including the high-p^ 
correction to the flow velocity due to the non-equilibrium 
effect. Next, the Raa and v-i of J/tp were calculated with 
this schematic model. Both the Cronin effect due to the 
initial transverse momentum broadening of gluons in nu- 
clei and the leakage effect due to the escape of initially 
produced J/ip from produced hot dense matter were in- 
cluded, and it was found that the Raa was suppressed at 
low pt and enhanced at high pt by both effects. How- 
ever, the leakage effect was much smaller than the Cronin 
effect. For the kinetic freeze-out temperature of J/tp, it 
was taken to be the same as the chemical freeze-out tem- 
perature, because the elastic cross section of J/ip is much 
smaller than that of light hadrons, which normally has 
a lower kinetic freeze-out temperature than the chem- 
ical freeze-out temperature. Finally, higher-order cor- 
rections in pQCD were included simply by multiplying 



FIG. 14: (color online) Elliptic flow V2 of J/ip in collisions at 
b = 9 fm with higher-order corrections. Solid, dashed, dot- 
ted, and dot-dashed lines correspond, respectively, to multi- 
plicative factors (A,B)=(1,1), (2,2), (5,4) and (10,6.5) for the 
higher-order corrections. The double-dot-dashed line is the 
V2 of J/tp if they are totally from regeneration. 



separate constant factors to the dissociation cross sec- 
tions of charmonia and to the elastic cross section of 
charm quarks. The former is related to the thermal de- 
cay widths of charmonia and the latter to the relaxation 
factor of charm quarks. Higher-order corrections were, 
however, not applied to the screening mass. We found 
that higher-order corrections enhanced both the dissoci- 
ation of charmonia and the relaxation of charm quarks. 
The former decreases the number of initially produced 
J/tp while the latter increases the number of regenerated 
J/ip. As a result, the total number of J/ip was found 
to change very little after including higher-order correc- 
tions. However, the fraction of initially produced J/ip 
and that of regenerated ones were affected significantly 
by higher-order corrections with larger higher-order cor- 
rections leading to a larger fraction of regenerated J/ip 
among total J /ip- This is the reason why the Raa of J/ip 
as a function of the number of participants was insensi- 
tive to the higher-order corrections, and both the two- 
component model and the statistical model can describe 
successfully the dependence of Raa on the number of par- 
ticipants. Applying same higher-order corrections to the 
transverse momentum dependence of the Raa of J/ip, 
we found that it was also insensitive to the fraction of 
regenerated J/ip among the total J/ip except at high px, 
as the small transverse momentum of regenerated J / ip is 
compensated by the large transverse flow velocity. On 
the other hand, the i>2 of J/ip is sensitive to the fraction 
of regenerated J / ip because the i>2 of initially produced 
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J /if) is almost zero. If the elastic cross section of J/ip in 
QGP is not small, the initially produced J/ip can then 
acquire nonnegligible v% through scattering with thermal 
partons. However, initially produced J/ip is expected to 
have a smaller binding energy and larger radius because 
of the high temperature. The loosely bound J/ip is ex- 
pected to be dissociated rather than scattered elastically 
when colliding with thermal partons. For regenerated 
J/ip, they were assume to be produced from thermalized 
charm quarks in QGP. It is not clear yet how close charm 
quarks are from thermal equilibrium in heavy-ion colli- 
sions. Recent experimental data on non-photonic single 
electrons, which are produced through the decay of heavy 
mesons, suggest that charm quarks are thermalized sig- 
nificantly [29(. If charm quarks are thermalized, result- 
ing thermally produced J/ip would share their thermal 
properties and thus acquire appreciable v%- A possible 
future analysis would be to include the effect of abrupt 
decrease in mass of the charmonia near T c (3l| . as de- 
creased mass will enhance the yield of regenerated J/ip 
and hence change its total contribution to V2- Presently, 
the V2 of J/ip has large error bars in the experimental 
data. If refined and precise data on the Vi of J/ip be- 
come available, they are expected to play an important 
role in discriminating between the two production mech- 
anisms for J/ip, i.e., initial production and regeneration. 
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Appendix 



Let the semimajor and the semiminor of an elliptic 
fireball produced in heavy- ion collisions be R y and R x , 
respectively, and the fluid velocities along these two axes 
be v y and v x , respectively, at a certain time during the 
evolution of the fireball. A short time At later, the 
semimajor becomes R' y = R y + v y At and the semimi- 
nor R' x = R x + v x At. Suppose that a certain point on 
the ellipse, (x/R x ) 2 + (y/R y ) 2 — 1, moves to a slightly 
expanded ellipse, (x/R x ) 2 + (y/R' y ) 2 = 1, in radial direc- 
tion keeping the same polar angle (p. Then the velocity 
of the point v{<p) satisfies 



-i 2 



{R(cp) +v(cp)At\ coscp 



R.i 



V r At 



! 2 



{R(<p) + v(<P) At} sin <P 



Ry + VyAt 



(34) 



where 



R(<f>) 



Rl 



Rl 



-1/2 



is the distance from the center of the ellipse to the point on the ellipse along the polar angle (p. Expanding Eq. (|3"4")l 
with respect to the infinitesimal time interval At gives 



{R(cP) cos^} 2 _ {R^smcf,} 2 + ^ 



Rl 



R 2 



v{<f>) 



2R(<P) cos 2 <P 2i?(</>) sin 2 <A 2R(cp) 2 v x cos 2 (p 2R 2 {cP)v y sin 2 <f>' 



Rl 



R 2 



Rl 



*5 



= 1, 



leading thus to the following expanding velocity of the 
surface of the fireball 



v(<t>) = R\<P) 



V x COS (p v y srn 



Rl 



(35) 



Because high-p^ particles are unlikely to follow the col- 
lective motion, both v and Av can be functions of pt- 
For simplicity, we modify only Av by multiplying it with 
the factor exp[— C(pt/ti)], where C is a fitting parameter 
301. 



Defining v = (v x +v y )/2 and Av = (v x —v y )/2, Eq. ([35 
can be rewritten as 



v(cP) =i? 3 (0) 



(v + Av) cos 2 <f> (v — Av) sin 2 



Rl 



R 3 



(36) 



15 



1] T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). [17 
2] M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92, 012001 

(2004). [18 
3] S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, Phys. 

Rev. D 69, 094507 (2004). [19 
4] R. L. Thews, M. Schroedter and J. Rafelski, Phys. Rev. 

C 63, 054905 (2001). [20 
5] M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. 

B410, 337 (1997). [21 
6] R. Arnaldi et al. [NA60 Collaboration], Nucl. Phys. A [22 

783, 261 (2007). 

7] A. Adare et al. [PHENIX Collaboration], Phys. Rev. [23 

Lett. 98, 232301 (2007). 
S] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. [24 

Lett. 96, 012304 (2006). 
9] B. Zhang, C. M. Ko, B. A. Li, Z. w. Lin and B. H. Sa, [25 

Phys. Rev. C 62, 054905 (2000). 
0] L. Grandchamp and R. Rapp, Nucl. Phys. A709, 415 [26 
(2002). 

1] L. Yan, P. Zhuang and N. Xu, Phys. Rev. Lett. 97, [27 
232301 (2006). 

2] A. Andronic, P. Braun-Munzinger, K. Redlich and [28 
J. Stachel, Nucl. Phys. A789, 334 (2007). [29 
3] T. Song, W. Park and S. H. Lee, Phys. Rev. C 81, 034914 

(2010). [30 
[14] X. Zhao and R. Rapp, Phys. Lett. B664, 253 (2008). 
[15] S. V. Akkelin, P. Braun-Munzinger and Yu. M. Sinyukov, [31 

Phys. Rev. C 81, 034912 (2010). 
[16] A. Adare et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 98, 232002 (2007). 



L. A. Linden Levy [PHENIX Collaboration], J. Phys. G 
36, 064013 (2009). 

Y. Park, K. I. Kim, T. Song, S. H. Lee and C. Y. Wong, 
Phys. Rev. C 76, 044907 (2007). 

T. Song, Y. Park, S. H. Lee and C. Y. Wong, Phys. Lett. 
B659, 621 (2008). 

M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95, 
122001 (2005). 

P. Levai and U. W. Heinz, Phys. Rev. C 57, 1879 (1998). 
P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. 
B465, 15 (1999). 

D. Teaney, J. Lauret and E. V. Shuryak, 
|arXiv:nucl th/0110037| 

R. J. Fries, B. Muller, C. Nonaka and S. A. Bass, Phys. 
Rev. C 68, 044902 (2003). 

S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 
69, 034909 (2004). 

S. Afanasiev et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 99, 052301 (2007). 

C. Silvestre [PHENIX Collaboration], J. Phys. G 35, 
104136 (2008). 

D. Krieg and M. Bleicher, Eur. Phys. J. A 39, 1 (2009). 
A. Adare et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 98, 172301 (2007). 

Y. Oh, Z. W. Lin and C. M. Ko, Phys. Rev. C 80, 064902 
(2009). 

K. Morita and S. H. Lee, Phys. Rev. Lett. 100, 022301 
(2008). 



